Effect of Two Level System Saturation on Charge Noise in Josephson Junction Qubits 
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We show that charge noise Sq in Josephson qubits can be produced by fluctuating two level 
systems (TLS) with electric dipole moments in the substrate using a flat density of states. At high 
frequencies the frequency and temperature dependence of the charge noise depends on the ratio 
J/J c of the electromagnetic flux J to the critical flux J c . It is not widely appreciated that TLS 
in small qubits can easily be strongly saturated with J/J c 2> 1. Our results are consistent with 
experimental conclusions that Sq ~ 1/f at low frequencies and Sq ~ / at high frequencies. 
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Noise and decoherence are a major obstacle to using 
superconducting Josephson junction qubits to construct 
quantum computers. Recent experiments 0, 0] indicate 
that a dominant source of decoherence is two level sys- 
tems (TLS) in the insulating barrier of the tunnel junc- 
tion as well as in the dielectric material used to fabricate 
the circuit. It is believed that these TLS fluctuators lead 
to low frequency 1/f charge noise Sq (/) d, 0, H, S 0] • 
However, at high frequencies, one experiment finds that 
the charge noise increases linearly with frequency Q. 
This has prompted some theorists to use a TLS den- 
sity of states linear in energy Q which is contrary to the 
constant density of states that has been so successful in 
explaining the low temperature properties of glasses such 
as the specific heat that is linear in temperature [loj . A 
linear distribution has been proposed in conjunction with 
a Cooper pair tunneling into a pair of electron traps ll| , 
and with electron hopping between Kondo-like traps to 
account for the charge noise 

However, these previous theoretical efforts have ne- 
glected the important issue of the saturation of the 
two level systems. Dielectric (ultrasonic) experiments 
on insulating glasses at low temperatures have found 
that when the electromagnetic (acoustic) energy flux 
J used to make the measurements exceeds the criti- 
cal flux J c , the dielectric (ultrasonic) power absorption 
by the TLS is saturated, and the attenuation decreases 
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explain the linear increase in the charge noise in Joseph- 
son junctions assumed that the TLS were not saturated, 
i.e., that J <c J c - This seems sensible since the charge 
noise experiments were done in the limit where the qubit 
absorbed only one photon. However, stray electric fields 
could saturate TLS in the dielectric substrate as the fol- 
lowing simple estimate shows. We can estimate the volt- 
age V across the capacitor associated with the substrate 
and ground plane beneath the Cooper pair box by setting 
CV 2 /2 — hui where Hu is the energy of the microwave 
photon. We estimate the capacitance C — e eA/L ~ 7 
aF using the area A = 40 x 800 nm 2 8] of the Cooper 
pair box, the thickness L = 400 nm of the substrate 



and the dielectric constant e = 10. Using u>/2ir = 10 
GHz, we obtain a voltage of V ~ 1.4 mV. The sub- 
strate thickness L of 400 nm yields an electric field of 
E ~ 3.4 X 10 3 V/m. For amorphous Si0 2 at / = 7.2 
GHz and Sir\L_ at / = 4.7 GHz, the critical rms voltage 
V c ~ 0.2 /j>V [2j], and with a capacitor thickness of 300 
nm, the critical field is E c ~ 0.7 V/m at T = 25 mK. So 
E/E c - 5 x 10 3 , and J/J c = (E/E c ) 2 - 2 x 10 7 > 1. 
A similar estimate shows that a single photon would 
even more strongly saturate resonant TLS in the insu- 
lating barrier of the tunnel junction, again resulting in 
J/J c 3> 1. However, there are only a few TLS in the 
oxide barrier of a small tunnel junction. (For a parallel 
plate capacitor with a specific capacitance of 60 fF //im 2 , 
L = 1.5 nm, A = 1 /im 2 and a dielectric constant of 
10, the volume is f2 = 1.5 x 10~ 21 m 3 . With a density 
of states P a ~ 10 45 / (Jm 3 ) ~ 663/ftGHz/mr 3 , there are 
only 2 TLS with an energy splitting less than 10 GHz.) A 
single fluctuator would have a Lorentian noise spectrum. 
The presence of 1// noise implies many more than 2 fluc- 
tuators. It is likely that these additional fluctuators are 
in the substrate. Our main point is that TLS in small 
devices are easily saturated. 

In this letter we explore the consequences of this sat- 
uration. We find that at high frequencies (Has 3> ksT), 
the frequency and temperature dependence of the charge 
noise depends on the ratio J/J c (oj, T) of the electromag- 
netic flux J to the critical value J c (to, T) which is a func- 
tion of frequency lu and temperature T. Starting from the 
fluctuation-dissipation theorem, we show that the charge 
noise is proportional to the dielectric loss tangent tan 5. 
We then calculate the dielectric loss tangent due to fluc- 
tuating TLS with electric dipole moments [3, 14, 17 1. At 
low frequencies we recover 1/f noise. At high frequen- 
cies tan S is proportional to 1/ yl + (J/Jjui, T)). In the 
saturation regime (J 3> J c (u>,T)), tani5 and hence, the 
charge noise, are proportional to y/ J c (ui,T)/ J. Some 
TLS experiments (lit [Til . [20| indicate that J c (u>,T) ~ 
u! 2 T 2 which implies that at high frequencies the charge 
noise and the dielectric loss tangent would increase lin- 
early in frequency if J 3> J c (w, T). Unlike previous theo- 
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retical efforts, we use the standard TLS density of states 
that is independent of energy, and can still obtain charge 
noise that increases linearly with frequency in agreement 
with the conclusions of Astafiev et al. Q. 

In applying the standard model of two level systems 
to Josephson junction devices, we consider a TLS that 
sits in the insulating substrate or in the tunnel barrier, 
and has an electric dipole moment p consisting of a pair 
of opposite charges separated by a distance d. The elec- 
trodes are located at z = and z = L and kept at the 
same potential. The angle between p and z-axis which 
lies perpendicular to the plane of the electrodes is 8 Q . 
The dipole flips and induces charge fluctuations on the 
electrodes. These induced charges are proportional to the 
z-component of the dipole moment, i.e., Q = \pcosdo/ L\. 

The TLS is in a double-well potential with a tunneling 
matrix element Ao and an asymmetry energy A [10| . The 
Hamiltonian of a TLS in an external ac field can be writ- 
ten as H = Hq + Hi, where Hq = ^{Aa z + Aq<7 x ), and 
H = — & z jy • ^ac{t)- Here cr xz are the Pauli spin matri- 
ces and £ac(i) = ^accosut is a small perturbing ac field of 
frequency u> that couples to the TLS electric dipole mo- 
ment. After diagonalization, the TLS Hamiltonian be- 
comes Hq = \E<j z , where E is the TLS energy splitting 
(i.e., E = \J A 2 + A 2 ) and, in the new basis, the perturb- 
ing Hamiltonian is Hi = —(a z A/E + o- x Aq/ E)p ■ £ a c(t)- 
The complete TLS Hamiltonian is similar to the Hamil- 
tonian Hs = — 7S • B for a spin 1/2 particle in a 
magnetic field given by B(t) = Bo + Bi(i), where the 
static field is — frjUo — (0,0, E) and the rotating field 
is fryBi(i) = {2A /E,0,2A/E)p • £ ac coswt [ID]. 7 is 
the gyromagnetic ratio and S = hcr/2. Therefore, the 
equations of motion of the expectation values of the spin 
components are given by the Bloch equations [22I ] with 
the longitudinal and the transverse relaxation times, Ti 
and T 2 . In the standard model, an excited two- level sys- 
tem decays to the ground state by emitting a phonon. 
The longitudinal relaxation rate, Tf" , is given by [10( 

E 



Tf 1 = aSA^coth 
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According to the Wiener-Khintchine theorem, the 
charge spectral density Sq(lo) is twice the Fourier trans- 
form ^q(ui) of the autocorrelation function of the fluc- 
tuations in the charge. From the fluctuation-dissipation 
theorem, the charge noise is given by: 



S Q (k,u>) = T 



4ft 
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where Q is the induced (bound) charge and 
Xq(fe, u>) is the Fourier transform of ^^(r : t]r' ,t') — 
{[Q(r,t),Q(r',t')]}/2h. We use Q = JP- dA, where P 
is the electric polarization density, and choose P z and 
dA\\z since Q ~ \p z \ to find ^{k.uo) = e A 2 xp z (k,uj), 
where £0 is the vacuum permittivity, A is the area 
of a plate of the parallel plate capacitor with ca- 
pacitance C, and x'p (fci w ) is the imaginary part of 
the electric susceptibility. Setting k = 0, and using 
EoXp (^) = £'(^) tan5(w), and C — e'A/L, we find 



Sq(u) = 



4HC 



-huj/k B T 



tan(5(w), 



(3) 



where Sq(uj) = Sg(k — 0,w)/fi, the volume of the ca- 
pacitor is f2 = AL, and s'(uj) = e' + etls{u) ~ e' = eoe r 
where e r is the relative permittivity. The frequency de- 
pendent £tls( w ) produced by TLS is negligible com- 
pared to the constant permittivity e' [To| | . 

The dynamic electric susceptibilities (x'(w), x"( w )); 
and hence the dielectric loss t ang ent, can be obtained by 

II [H H3. (Shnirman et 



where a = [jj/ fapft 4 )] [(l/cf ) + (2/cf)] where p is the 
mass density, Cg is the longitudinal speed of sound, c t 
is the transverse speed of sound, and 7^ is the defor- 
mation potential. The distribution of TLS parameters 
can be expressed in terms of E and T-y: P{E,T\) = 
P /(2Tiy / l-T mm {E)/Ti) [13, [H where P is a con- 
stant. The minimum relaxation time T m i n (E) corre- 
sponds to a symmetric double-well potential (i.e., E — 
Ao). The transverse relaxation time Ti represents the 
broadening of levels due to their mutual interaction 24 1 . 

General Expression for Charge Noise: We derive a 
general expression valid at all frequencies for the charge 

noise in terms of the dielectric loss tangent tan^(cj) = conducting qubits is often extremely short, t 
e" (lu) I e' (uj) , where e'(ui) and e"{tu) are the real and imag- 
inary parts of the permittivity. 



solving the Bloch equations 
al. [9(] gave a reformulation of the Bloch equations for 
TLS in terms of density matrices and the Bloch- Redfield 
theory.) The electromagnetic dispersion and attenuation 
due to TLS has two contributions. First there is the re- 
laxation process due to the modulation of the TLS energy 
splitting by the incident wave resulting in readjustment 
of the equilibrium level population. This is described by 
Xz{u) that comes from solving the Bloch equations for 
S z which is associated with the population difference of 
the two levels. The second process is the resonant ab- 
sorption by TLS of phonons with Hlo = E. Resonance 
is associated with x±( UJ ) = Xx(w) ± ix y {to) since a x and 
<j y are associated with transitions between the two lev- 
els. The total dielectric loss tangent is the sum of these 
two contributions: tan(5 = tancS^i + tan 5res- The 
steady-state solution of the Bloch equations and the re- 
sulting dielectric loss tangent for TLS are well known 

EE El, E| . 

In the steady-state regime the experimental values for 
the relaxation times Ti and T2 are considered small com- 
pared to the electromagnetic pulse duration, t p . One 
might ask if one should use the transient solution 21 1 
of the Bloch equations since the pulse applied to super- 
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sec Q. We find that the transient z-component of 
the magnetization [2l[, S z (t), decays exponentially to 
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the equilibrium value denoted by S z eq , i.e.: S%(t) — 
S° eq + exp(-t/T*)[S°(0)-S° ieq ], where S° z (0) is the ini- 
tial value of S" and S° zeq = -htanh(E/2k B T)/2. The 
transient relaxation time T* is given by 



1 + (J/ J c (w, T)) x g(ci>, lli q , T 2 ) ' 
g(w,w ,T 2 ) 



where 
1 



l + T 2 (w-w ) 2 l + T|(w + a;o) 2 



Here J/J c (u,T) = (p' -^JhfT^ p' = (A /e)p repre- 
sents the induced TLS dipole moment, and lvq = "/B = 
—E/h. In the saturated regime, for J 3> J c (u>,T), we 
find that T* w [T 2 (p' • W*) 2 ] -1 - Using p' = 3.7 D, 
= 3.4 x 10 3 V/m, and T 2 = 8 at T = 0.1 K 
(l9L I2II ] yields T* ~ 8 x 10 -13 sec at resonance when 
u/2tt ~ luo/2tt — 10 GHz. This value is much shorter 
than the typical pulse length used in the Josephson junc- 
tion qubits experiments and therefore the results from 
steady-state saturation theory can be used. However, 
in the unsaturated regime, T* w T\ » 8 x 10" 8 sec at 
T = 0.1 K, so transient effects can be important. 

High Frequency Charge Noise: At high frequencies 
(HF) (Hu> ^> k B T) the dielectric loss tangent is dom- 
inated by resonant (RES) absorption processes (i.e., 
tanS HF -tan 5 RES ) EEQUl: 



tan 5hf(w, T) 



3e' 



tanh(fiw/2fc B T) 



1 



j/Jc{u,Ty 

(5) 

where J c (cj,T) = 3h 2 e'v/(2p 2 T 1 T 2 ) and t> is the speed 
of light in the solid. Eq. ([5]) comes from integrating 
over the TLS distribution [14| . However, if no inte- 
gration is done due to the small number N n of TLS, 
e.g., in the tunnel junction barrier, tan S h f (oj , T) — 
N p 2 T 2 /(3s'Ml) xtanh(hw/2k B T)[\ + J/J c (w, T)]" 1 . So 
for high intensities the (J/Jc)" 1 / 2 dependence of tan<5 
becomes a (J/J c ) _1 dependence. The frequency and 
temperature dependence of tan 5, and hence of Sq(u>), 
depends on J/ J c (u>,T). 

At low intensities (J <C J c (lo,T)) in the unsaturated 
steady-state resonant absorption regime, the dielectric 
loss tangent is constant: 



tan 5 



HF 



tt P 2 P / (3s') 



(6) 



For e r = 10, P « 10 45 (Jm 3 )- 1 [27|, |28|, and p = 3.7 
D (which corresponds to the dipole moment of an OH" 
impurity [111), we estimate tanSnF ~ 1.8 x 10" 3 . This 
result agrees well with the value of 6 ~ 1.6 x 10" 3 re- 
ported in Ref. In this regime the charge noise is 
constant: Sq/e 2 ~ AirhCp 2 P / '(3eV). For C = 7 aF, 
S Q /e 2 ~2x 10" 16 Hz" 1 . 

For high field intensities J ^> J c (uj,T), we obtain 
tan^F = 7rp 2 P /(3e') x ^/J c (w,T)/J. The J" 1 / 2 de- 
pendence of tan 5 has been found for materials such as 



amorphous Si0 2 and amorphous SiN^ Q]. In the satu- 
rated resonant absorption regime the charge noise is given 
by 



S Q (u,T) 4hCTT P 2 P J c (w,T) 



3e' 



J 



(7) 



So the frequency and temperature dependence of the 

noise is determined by Z\ and T 2 : Sq(uj,T) ~ 
w 2 

(TiT 2 ) . Experiments and theory 
~ T m where m ranges from 1 to 2.2 
31 , H, 33] . T 2 decreases with increas- 



(4) y/J c (u,T) 



find that T,~ 
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ing frequency 2fJ but the exact frequency dependence 
is not known, so we will ignore it in what follows. T\ 
in the symmetric case with hw = E = A D is given by 



Eq. ©: T{ 



UJ 



implying that at high frequencies 



and intensities Sq(uj : T) ~ w 3 / 2 T m / 2 , where m/2 varies 
between 0.5 and 1.1. If there are only a very few TLS, 
then Sq(w,T) ~ J c (u,T) ~ (TiT 2 ) - u 3 T m . Although 
the experimental frequency dependence was reported to 
be linear Q, the scatter in the data is large enough to 
allow for a steeper frequency dependence. In fact, the 
dependence is much steeper for data at the degeneracy 
point. 

Experimental measurements [l5| on Si0 2 at / = 10 
GHz and 0.4 < T <1 K find J c (w,T) = 25 mW/cm 2 x 
(T/0.4K) 4 , and that J/J c (u,T) varies between 10" 2 and 
10 4 . This implies that the charge noise Sq/e 2 should 
vary between 2 x 10" 16 Hz" 1 and 2 x 10" 18 Hz" 1 at 
T = 0.1 K. However, other measurements have found 
JAuj,T) ~ to n T 2 [3, where n is equal to either 0, 
Q or 2 [HIU. If J c (w,T) ~ w 2 , then tan5(u) - w, 
and Sq(lo) ~ lu at high frequencies and high intensities, 
which agrees with the recent experiments by Astafiev et 
al. [8j . It would be interesting measure the temperature 
dependence of Sq experimentally. 

While currently there are no direct experimental mea- 
surements of the HF charge noise, in Ref. [8j Sq(oj) has 
been deduced by measuring the qubit relaxation rate Ti 
versus the gate induced charge q for a Cooper pair box 
with a capacitance Cb, Josephson energy Ej (in the GHz 
range), and electrostatic energy U — 2eq/Cb- From Ti = 
■KS u (iv)sm 2 {9)/2h 2 , where sin0 = E 2 /(E 2 +U 2 ) and hu) 

1 /2 

equals the qubit energy splitting (hu = (U 2 + E 2 ) ), 
and Su(u) — (2e/Cb) 2 S q (uj), they obtained the charge 
noise S q (uj) at high frequency [1]. We can compare our 
results with these experiments by reversing this proce- 
dure to find Ti(q) from S q (uj). We find for saturated TLS 
that Ti at the maximum (q = 0) is of order 10 8 s _1 [35| 
and increases as E 2 in good quantitative agreement with 
the experimental results from Fig. 2 and 3 of Astafiev et 
al. [|. We find that I\ 

max: the maximum value of Ti 
(at q — 0), increases with the frequency / = Ej /h in the 
saturated regime but is independent of frequency at low 
intensities J <C J c (w, T). 

Low Frequency Charge Noise: We now show that we 
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can recover the low frequency 1/f charge noise using 
Eq. ([3]). At low frequencies (LF) where Huj <C k B T, only 
the relaxation absorption process contributes: tanS^p ~ 
tan<W = np 2 P /(6e') O. Eq. © gives [3il 



SqU) 2k B T 1 1 T (P\ 2 1 



2 /2C 



2tt/ 



(8) 

To estimate the value of <Sq, we use p = 3.7 D, 
P » 10 45 (Jm 3 )' 1 , L = 400 nm, and A = 40 x 800 
nm 2 . At T = 100 mK and / = 1 Hz, we obtain 
Sq/c 2 = 2 x 10~ 7 Hz -1 , which is comparable to the ex- 
perimental value of 4 x 10 -6 Hz -1 deduced from current 
noise [?]|. As Eq. ([5]) shows, the standard TLS distri- 
bution gives low frequency 1/f charge noise that is lin- 
ear in temperature, while experiments find a quadratic 
temperature dependence 0, 0] • This implies that at low 
frequencies and temperatures contributions from other 
mechanisms may dominate the charge noise 



11 



To conclude, we have shown that the frequency and 
temperature dependence of high frequency charge noise 
in Josephson junction devices depends on the ratio 
J/J c (u>, T) of the electromagnetic flux to the critical flux. 
Using the standard theory of two level systems with a 
flat density of states, we find that the charge noise at 
high frequencies can increase linearly with frequency and 
temperature if J/J c (u,T) ;§> 1. This agrees with the 
conclusions of recent experiments on the high frequency 
charge noise in Josephson junction qubits [8[ which our 
estimates show are in the strongly saturated limit. 
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